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FOREWORD 


This  report  presents  results  from  an  analysis  of  Global  Positioning  System  (GPS)  satellite 
and  marine  acoustic  observations  to  extend  WGS  84  geodetic  control  onto  the  seafloor.  The  data 
were  collected  as  part  of  an  experiment  in  marine  geodynamics  conducted  by  the  United  States 
Geological  Survey  (USGS)  with  support  from  the  University  of  Hawaii;  the  Naval  Postgraduate 
School;  the  Defense  Mapping  Agency  (DMA);  the  Naval  Surface  Warfare  Center,  Dahlgren 
Division  (NSWCDD);  and  the  U.S.  Navy  Submarine  Group  One.  At  NSWCDD,  data  analysis 
was  performed  in  the  Space  and  Surface  Systems  Division  of  the  Strategic  and  Space  Systems 
Department  imder  sponsorship  of  the  Defense  Mapping  Agency. 

TTie  authors  would  like  to  acknowledge  Dr.  Janet  Morton  of  the  USGS  in  Reston,  Virgima, 
for  her  dedication  in  organizing  the  experiment;  Dr.  Muneendra  Kumar  of  the  DMA  who 
recognized  the  geodetic  opportunity  represented  by  the  experiment;  Ms.  Carol  Reiss  of  the  USGS 
in  Menlo  Park,  California,  for  providing  GPS  and  acoustic  data  sets  and  for  explaining  maiiy 
details  associated  with  the  USGS  data  acquisition  and  processing;  Professor  Steven  Tucker  of  the 
Naval  Postgraduate  School  for  providing  the  sound  velocity  profile  used  for  refraction  correction; 
and  Mr.  Paul  Jackins  of  NSWCDD  who  provided  software  for  acoustic  ray  tracing. 

This  report  has  been  reviewed  by  J.  L.  Sloop,  Head,  Space  and  Surface  Systems  Division. 


Approved  by: 

R.  L.  SCHMIDT,  Head 

Strategic  and  Space  Systems  Department 
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ABSTRACT 


The  results  of  an  analysis  of  a  complex  data  set  acquired  during  the  United  States 
Geological  Survey’s  Marine  Crustal  Deformation  Study  are  presented.  The  experiment,  which 
commenced  in  the  spring  of  1992  in  a  region  of  the  Pacific  known  as  the  Juan  de  Fuca  ridge, 
represents  a  first  attempt  to  locally  monitor  plate  dynamics  in  the  marine  environment  using  a 
network  of  tripod-mounted,  dual-frequency  acoustic  transponders.  The  aim  is  to  collect  over  a 
period  of  years  time-series  measurements  of  extension  rates  along  the  southern  Juan  de  Fuca  ridge. 
In  addition,  by  collecting  a  combination  of  Global  Positioning  System  satellite  tracking  data,  low- 
frequency  acoustics  data,  and  water  column  pressure,  conductivity,  and  temperature  at  depth,  it 
was  possible  to  extend  geodetic  control  from  land  onto  the  seafloor.  The  methods  to  accomplish 
this  latter  goal  of  the  experiment  are  described,  as  well  as  the  final  results.  Baseline  comparisons 
between  several  of  the  solutions  obtained  during  this  analysis  are  presented  along  with 
recommendations  for  additional  data  collection  and  analysis. 
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INTRODUCTION 


The  long-term  goal  of  extending  geodetic  control  into  the  oceans  is  moving  closer  to 
becoming  a  practical  reality  due  to  improvements  in  acoustical  instrumentation  and  through  the  use 
of  the  Globd  Positioning  System  (GPS).  Recent  experimentation  in  199 1  by  the  Jet  Propulsion 
Laboratory  and  the  Scripps  Institution  of  Oceanography,  University  of  California  at  San  Clego, 
followed  in  1992  by  an  experiment  organized  by  the  United  States  Geological  Survey  (USGS), 
demonstrates  that,  although  not  by  any  means  routine,  these  combined  measurement  technologies 
enable  geodetic  positioning  under  the  seas  to  be  obtained  in  a  global  reference  frame. 

During  the  month  of  May  1992,  an  acoustic  transponder  network  was  deployed  across  the 
Juan  de  Fuca  ridge,  300  nmi  off  the  Oregon  coast.  This  experiment  sponsored  by  the  USGS  and 
the  Defense  Mapping  Agency  (DMA)  has  a  principal  goal  of  monitoring  the  crustal  spreading  rate 
in  the  marine  environment  through  long-term  observation  at  a  fixed  seafloor  network.  A  second 
objective  of  the  experiment  was  to  determine  geodetic  positions  for  the  deep  ocean  acoustic 
transponder  network  in  a  global  reference  frame,  thus  extending  WGS  84^  geodetic  control  onto 
the  seafloor.  The  experiment  represents  a  partnership  that  includes  contributions  from  the  USGS, 
the  University  of  Hawaii;  the  Naval  Postgraduate  School  (NPS);  the  DMA;  the  Naval  Surface 
Warfare  Center,  Dahlgren  Division  (NSWCDD);  and  the  U.S.  Navy  Submarine  Group  One. 


GEODETIC  POSITIONING  IN  THE  MARINE  ENVIRONMENT- 
EXPERIMENT  PLANNING 


The  main  focus  of  the  USGS  experiment  is  to  acoustically  measure,  over  time,  crustal 
dynamics  under  the  sea.  However,  NSWCDD  and  DMA  participated  with  the  USGS  in  the  design 
of  this  experiment  because  of  an  interest  in  estimating  geodetic  coordinates  for  seafloor 
instrumentation  using  a  combination  of  low-frequency  acoustics  and  GPS  satellite  measurements. 

Prior  to  going  to  sea,  an  analysis  was  carried  out  to  investigate  the  potential  benefits  of 
combining  GPS  satellite-based  positioning  and  acoustics  data  for  determining  the  geodetic 
coordinates  of  points  on  the  ocean  floor.  This  simulation  study  provided  indications  of  the 
accuracy  that  might  be  expected  for  seafloor  positioning,  as  well  as  a  design  for  shipboard  data 
acquisition.  Principally,  this  analysis  was  peif^ormed  using  the  Ohio  State  University  Geometric 
and  Orbital  Program^  software  developed  at  the  Ohio  State  University’s  Department  of  Geodetic 
Science  and  Surveying,  modified  by  NSWCDD  for  this  marine  application. 

This  ship-track  covariance  analysis  was  undertaken  to  determine  how  to  maximize  the 
geometric  strength  (figure  of  merit)  of  the  low-frequency  acoustic  data  acquisition.  Four 
transponders  were  considered  in  this  simulation  in  the  approximate  locations  of  the  planned  USGS 
deployment  located  at  the  comers  of  a  square  1  km  on  a  side  at  a  22()0-m  depth.  Three  ship-track 
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patterns  were  analyzed  in  determining  an  acoustics  data  acquisition  strategy.  Each  pattern 
consisted  of  a  square  (see  Figure  1)  augmented  with  north-south  and  east- west  crossing  tracks 
lying  within  its  boundary,  centered  within  the  square.  The  side  lengths  for  the  concentric  squares 
were  1,  3,  and  5  km.  The  squares  were  centered  over  the  transponder  array.  In  each  case,  the 
simulated  acoustic  range  data  acquired  along  the  track  pattern  (square  plus  crossing  tracks) 
provided  135  equally  spaced  events.  An  event  is  defined  as  a  set  of  four  acoustic  range 
measurements  between  the  ship  and  the  transponder  network  obtained  while  GPS  data  are  collected 
to  at  least  four  satellites,  completing  (conceptually)  a  double  pyramid  (see  Figure  2).  Thus  data 
density  was  spacially  greatest  for  the  1-km  square  pattern  and  least  for  the  5-km  square.  The 
designations  for  the  test  patterns  were  pattern  1  for  the  3-km  square,  pattern  2  for  the  5-km  square, 
and  pattern  3  for  the  1-km  inner  square.  A  data  simulator  was  developed  that  generated  acoustic 
ranges  to  the  transponders  from  a  ship  that  traversed  each  track  pattern. 


HGURE  1.  SHIP  TRACK  PATTERNS  FOR 
LOW-FREQUENCY  DATA  ACQUISITION 


These  simulated  acoustic  observations  were  then  processed  using  a  least-squares  estimation 
technique  to  determine  the  affect  of  ship-track  geometry  on  position  uncertainty  for  the  array. 
Results  were  obtained  first  for  pattern  1  using  different  constraints  on  the  accuracy  of  the 
associated  ship  positions  that  would  be  determined  at  sea  using  GPS  satellite  measurements.  The 
accuracy  to  which  a  ship  can  be  positioned,  both  absolute  and  relative,  varies  greatly  depending  on 
GPS  data  types,  collection  sites,  and  processing  modes.  Accuracies  of  10  m,  1  m,  and  10  cm  were 
considered  in  this  study.  Acoustic  ranges  were  assumed  to  have  a  1-m  standard  error  for  those 
slant  ranges  presented  by  the  ship-to-array  geometry.  Figure  3  provides  the  latitude,  longitude, 
and  ellipsoidal  height  uncertainty  (1  sigma)  for  transponder  T1  as  a  function  of  ship  position 
uncertainty.  Results  for  transponders  T2,  T3,  and  T4  were  almost  identical.  These  results 
indicated  that  ship  positioning  accuracy  of  1  m  relative  to  terrestrial  geodetic  control  would  be  more 
than  sufficient  for  the  experiment.  Improving  position  accuracy  using  more  involved  GPS  data 
acquisition  and  processing  would  improve  transponder  positioning  only  marginally,  this  limitation 
being  due  to  the  accuracy  of  the  low-frequency  acoustic  ranges. 
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sea  surface 


Subsequent  simulations,  each  with  135  events,  were  performed  for  track  patterns  2  and  3 
however  using  only  1-m  GPS  ship  position  uncertainty.  Figure  4  provides  these  additional 
results.  Pattern  2  (case  2)  provided  improved  geometric  strength  for  array  latitude  and  longitude 
over  track  pattern  1,  but  height  uncertainty  remained  about  the  same.  Acoustic  data  simulated 
along  pattern  3  (case  3)  provided  a  weaker  solution  than  either  patterns  1  or  2  with  height 
uncertainty  affected  least.  Finally,  Figure  4  provides  the  results  when  data  from  all  three  track 
patterns  were  combined  (case  4).  These  405  events  in  the  combined  solution  (1620  acoustic 
ranges)  were  distributed  with  a  higher  density  close  to  the  array.  Observation  density  then 
decreased  as  the  distance  from  the  array  increased.  The  combined  result  gave  standard  errors  for 
the  array  of  approximately  16,  28,  and  10  cm,  respectively,  for  the  latitude,  longitude,  and 
ellipsoidal  height  of  each  transponder  instrument.  The  difference  in  the  latitude  and  longitude 
uncertainties  was  a  result  of  placing  the  GPS  position  uncertainty  of  1  m  on  the  Cartesian 
coordinates  of  the  ship’s  positions  rather  than  on  geodetic  coordinates.  The  combined  result  is 
better  than  transponder  uncertainty  resulting  from  the  use  of  10-cm  GPS  positioning  with  data 
from  any  one  single  track  pattern  having  35  events  (compare  Figures  3  and  4).  Thus  an  increase  in 
acoustic  data  acquisition  using  multiple  tracks  seemed  preferable  to  obtaining  the  most  accurate 
GPS  positioning  possible  and  would  reduce  requirements  for  more  complex  GPS  data  acquisition 
and  processing. 
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10  1  0.1 
GPS  POSTION  UNCERTAINTY  (M) 


HGURE  3.  TRANSPONDER  ARRAY  UNCERTAINTY  VS.  GPS 
POSITIONING  UNCERTAINTY  FOR  TRACK  PATTERN  1 


RGURE  4.  TRANSPONDER  UNCERTAINTY  ASSOCIATED  WITH  1-M  GPS 
SHIP  POSITIONING;  CASE  (I):  PATTERN  I,  CASE  (2);  PATTERN  2, 
CASE  (3):  PATTERN  3,  CASE  (4);  COMBINED  SOLUTION 
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Track  patterns  1,  2,  and  3  together  represent  a  total  of  46  km  of  ship  track.  With  some 
redundancy  and  extra  tracking  to  connect  the  three  patterns,  the  entire  track  set  could  be  collected  in 
6  hr  if  the  ship  were  to  transit  the  approximately  53-km  distance  at  an  average  speed  of  8.8  km/hr. 
However,  restrictions  on  the  acoustic  data  acquisition  speed  of  the  ship  would  need  to  be  taken  into 
account  for  fmal  planning. 

GPS  availability  during  the  data  collection  period  would  be  able  to  support  continuous  four 
or  four-plus  satellite  coverage  for  relative  positioning  of  the  ship.  Ship  position  accuracy  as  a 
function  of  time  for  the  available  GPS  satellite  coverage  however  was  not  considered  in  this 
analysis.  It  was  believed  that  sufficient  GPS  coverage  would  be  available  during  the  data 
acquisition  phase  of  the  experiment  to  obtain  sufficient  accuracy  relative  to  a  fixed  land  site  at 
Newport,  Oregon.  There  were  no  requirements  to  collect  all  acoustic  and  GPS  data  in  a  single 
session.  TTie  ship  collection  could  be  performed  in  segments  over  several  days. 

These  study  results  indicate  the  potential  power  of  combining  GPS  observables  and 
acoustic  ranges  to  obtain  undersea  geodetic  positions  for  points  of  interest  It  is  believed  that  this 
approach  represents  the  most  likely  future  technology  to  extend  geodetic  control  into  the  ocean 
environment.  However,  this  study  did  not  consider  the  effect  of  systematic  errors  such  as  in  GPS 
ship  positioning  or  in  acoustic  refraction  correction.  It  can  therefore  only  provide  (optimistic) 
results  reflecting  the  geometric  contributions  to  such  a  solution  as  a  guide  for  data  acquisition. 


uses  MARINE  GEODYNAMICS  EXPERIMENT 


EXPERIMENT  BACKGROUND 

The  USGS  deployed  four  dual-frequency  Oceano  SI-43 1-3  acoustic  transponders  in 
tripods  southwest  of  an  area  containing  low-frequency  Benthos  transponders  previously  deployed 
by  the  National  Oceanic  and  Atmospheric  Administration  (NOAA).  The  site  selected  for  the 
deployment  was  300  nmi  west  of  Astoria,  Oregon  at  a  point  along  the  Juan  de  Fuca  ridge  where 
crustal  spreading  rates  were  estimated  to  be  6  cm  per  year.  Of  the  four  transponders  deployed, 
three  functioned  properly,  providing  an  extensive  acoustics  data  set  over  a  two-day  calibration 
period.  Acoustics  data  were  also  collected  from  several  of  the  NOAA  instruments.  Figure  5 
provides  the  locations  of  the  four  USGS  instruments  and  the  location  of  one  NOAA  transponder 
whose  data  were  analyzed  in  this  study.  Notice  in  Figure  5,  and  in  the  following,  that  the 
transponder  designations  T5,  T6,  T7,  and  T4  are  consistent  with  USGS  data  files  and  not  with 
notations  used  in  the  previous  section  on  experimental  planning. 

The  experiment,  as  it  pertained  to  precise  seafloor  positioning,  consisted  of  collecting 
several  types  of  data  over  a  common  time  interval,  then  using  this  data  to  transfer  geodetic  control 
from  land  onto  the  seafloor.  Figure  6  provides  the  general  layout  of  the  experiment  data  collection. 
Details  on  the  specific  data  sets  follow. 
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HGURE  6.  DATA  ACQUISITION  DURING  USGS  JUAN  DE  FUCA  RIDGE  EXPERIMENT 
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RELATIONSHIP  TO  EXPERIMENT  PLANNING 

Although  a  large  data  set  was  acquired  during  the  calibration  period,  after  collation  and 
editing,  the  actual  geometry  of  the  remaining  ship  tracks  above  the  USGS  array  covered  an  area 
most  analogous  to  track  pattern  3  considered  earlier.  Therefore,  the  data  set  used  in  the  latter  parts 
of  this  study  lacked  the  geometrical  strength  associated  with  the  broader  collection  patterns 
previously  studied.  This  issue  will  be  considered  when  the  results  from  the  data  analysis  are 
interpreted. 


INSTRUMENTATION  AND  DATA  COLLECTION 


Data  collected  aboard  the  Navy  Deep  Submergence  Vehicle  Turtle  and  Support  Ship  Laney 
Chouest  during  the  site  visit  included  dual-  and  single-frequency  pseudorange  and  carrier  phase 
observations  from  shipboard  Ashtech  LDXII  and  3DF  GPS  receivers  providing  ship  location  and 
three-dimensional  ship  attitude,  low  frequency  (8-16  kHz  band)  acoustics  between  the  vessel  and 
the  transponder  networic  providing  round-trip  travel  times  between  the  survey  ship  and  the  USGS 
and  NOAA  transponders,  and  direct  measurements  of  electrical  conductivity,  temperature,  and 
pressure  within  the  water  column  to  determine  a  sound  velocity  profile.  Because  of  a  failure  of  the 
Oceano  RM-201  range  meter,  high  frequency  (50-108  kHz  band)  intersite  acoustic  range  data 
acquisition  was  commanded  between  the  USGS  transponders;  however,  data  retrieval  was  not 
possible.  This  latter  data  will  eventually  be  retrieved  to  provide  direct  measurement  of  crustal 
dynamics.  It  will  also  serve  to  validate  the  geodetic  positioning  results  from  this  analysis.  In 
addition,  a  fixed  GPS  land  site,  station  VENTS,  at  the  Hatfield  Marine  Science  Center  in  Newport, 
Oregon,  was  occupied  with  a  dual-frequency  Ashtech  LDXn  GPS  receiver  to  support  differential 
positioning  of  the  ship  relative  to  known  geodetic  control.  Although  data  were  collected  over 
several  days,  the  data  analysis  was  confined  to  calibration  lines  that  spanned  a  period  that  began  at 
hour  18  on  day  142  and  finished  at  the  end  of  hour  3  on  day  143.  Because  the  different  data  types 
were  not  consistently  collected  over  the  same  periods,  the  entire  acoustics  data  set  could  not  be 
fully  exploited.  Figure  7  illustrates  the  time  intervals  during  days  142  and  143  when  various  data 
sets  were  collected. 


nOURE  7.  CALBRATION  TRACKS  FOR  JUAN  DE  FUCA  RIDGE  EXPERIMENT- 
DATA  COLLECTION  INTERVALS  ON  DAYS  142  AND  143 
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PRELIMINARY  DATA  PROCESSING  AND  DATA  EDITING 


After  receipt  of  the  data  sets  collected  aboard  the  Laney  Chouest  and  the  GPS  tracking  data 
acquired  at  Newport,  Oregon,  several  data  processing  steps  were  undertaken.  These  included 
editing  acoustics  data  to  remove  multiple  returns,  processing  GPS  receiver  range  measurements  to 
position  the  ship  relative  to  known  geodetic  control  on  land,  examining  acoustics  data  residuals  to 
eliminate  questionable  data,  evaluating  interpolating  schemes  to  determine  ship  location  at  acoustic 
signal  receive  times,  connecting  GPS  antenna  locations  to  the  transducer  through  the  use  of  GPS- 
derived  ship  attitude  data,  and  editing  data  sets  to  produce  complete  data  coverage  over  selected 
time  interv^s.  Several  of  these  procedures  were  executed  at  the  beginning  of  the  data  processing, 
others  were  added  as  the  data  analysis  proceeded,  refining  methods  and  including  more  complete 
modeling.  The  following  paragraphs  summarize  these  procedures. 


INITIAL  EDITING  TO  REMOVE  MULTIPLE  RETURNS 

Each  transponder  within  the  USGS  and  NOAA  arrays  operates  with  a  unique  response 
frequency.  The  frequencies  of  the  transponders  whose  data  were  processed  in  this  study  are  given 
in  Table  1,  along  with  their  approximate  positions.  During  a  first  evaluation  of  the  acoustics  data 
collected  aboard  ship,  it  was  readily  evident  that  the  data  files  consistently  contained  multiple  (up  to 
four)  returns  from  each  transponder.  This  was  clearly  a  result  of  multiple  reflections  of  the 
transmitted  acoustic  signal  off  the  seafloor  and  ocean  surface,  which  then  triggered  additional 
chirps  from  the  transponder  instrumentation. 

Using  approximate  locations  for  the  transponder  and  ship  positions,  an  initial  editing  of  all 
acoustic  data  was  completed  to  identify  and  eliminate  all  but  principal  returns  from  each 
transponder  associated  with  a  20-sec  transducer  pulse  rate.  These  edited  files,  configured  on  an 
hourly  basis,  represented  the  source  data  for  all  subsequent  processing  steps. 


TABLE  1.  TRANSPONDER  INFORMATION 


Instrument  Number 

Latitude* 

(deg) 

Longitude* 

(dts) 

NOAA  04 

44.706 

229.679 

-1987 

13000 

USGS  05 

44.667 

229.635 

-2205 

11500 

USGS  06 

44.671 

229.647 

-2205 

12000 

USGS  07 

44.673 

229.639 

-2205 

12500 

"‘Apinuxiinate  position  in  WGS  84 
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GPS  SfflP  POSITIONING 

GPS  ship  positioning  was  accomplished  using  the  Ashtech  software  program  PPDIFF, 
which  was  designed  to  differentially  post-process  navigation  solutions  collected  in  the  field.  Since 
an  initially  known  baseline  was  not  established,  an  alternate  processing  program  GPPS  for  the 
kinematic  mode  was  not  used.  GPPS  would  have  been  more  desirable  from  an  accuracy 
standpoint  because  it  could  have  used  precise  GPS  ephemerides.  PPDIFF  required  the  raw 
measurements  that  were  downloaded  from  the  GPS  receivers  located  at  the  base  station  and 
onboard  the  ship,  as  well  as  the  broadcast  GPS  ephemerides. 

The  accuracy  of  the  results  from  PPDIFF  were  expected  to  range  from  5  to  15  m  under 
optimal  conditions  over  the  distances  present  in  this  experiment.  Errors  due  to  positioning  were 
expected  to  be  systematic  (correlated)  because  of  errors  in  the  GPS  ephemerides. 


RESIDUAL  EDITING 

Acoustic  data  sets  known  as  calibration  lines  were  placed  in  single  hour  files  beginning  on 
day  142,  hour  18  and  ending  after  day  143,  hour  3.  The  data  files  initially  contained  a  ship 
position  for  each  signal  transducer  pulse  time  (20-sec  epochs)  in  a  local  reference  frame.  Ship 
positions  in  this  frame  were  estimated  using  acoustics  data  collected  from  four  of  the  NOAA 
transponders  whose  coordinates  in  this  local  frame  were  previously  defined.  The  local  frame  was 
defined  by  the  directions  east,  north,  and  down,  with  its  origin  on  the  ocean  surface. 

During  data  processing,  it  was  necessary  to  edit  out  spurious  data  points  and  to  determine 
appropriate  weighting  for  the  data.  This  was  done  by  performing  initial  adjustments  of  ship  and 
transponder  coordinates  and  computing  means  and  standard  deviations  for  the  residuals  of  fit. 
Once  computed,  these  statistics  served  as  a  basis  for  refining  the  weighting  of  the  data  and  as  a 
means  to  edit  (remove)  data  values  whose  residuals  exceeded  a  linear  combination  of  the  mean  and 
a  selected  multiplier  of  the  standard  deviation.  Since  a  systematic,  unmodeled  signal  was  generally 
present  in  the  residuals  of  each  hourly  data  set,  the  multiplier  used  to  edit  data  in  each  adjustment 
was  heueristically  determined.  Examples  of  these  editing  criteria  are  presented  in  the  following 
paragraphs. 


SHIP  POSITION  INTERPOLATION 

Acoustic  data  were  collected  using  a  20-sec  transducer  pulse  rate.  The  principal  returns 
from  the  transponders  occurred  within  3  to  6  sec  after  the  transducer  pulse  event  depending  on  ship 
location.  After  completion  of  GPS  data  processing  to  determine  ship  position,  acoustic  data  files 
were  modified  to  include  interpolated  GPS  ship  positions  at  pulse  emit  times.  These  GPS 
positions  in  WGS  84  were  determined  using  an  eight-point  Lagrangian  interpolation  routine  that 
operated  on  the  GPS  ship  position  files.  The  interpolated  values  were  then  substituted  into  the 
acoustic  data  files  for  positional  information  in  the  local  reference  frame  as  provided  by  the  USGS. 
These  GPS  ship  positions  were  thus  developed  at  the  same  resolution  as  the  original  acoustics 
data,  not  at  the  1-sec  resolution  of  the  actual  GPS  data. 
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Although  this  reduced  the  number  and  size  of  the  working  data  files,  several  disadvantages 
resulted  from  this  procedure  or  were  present  due  to  other  circumstances.  First,  the  schedule  of 
operation  of  the  dual-frequency  GPS  receiver  on  the  ship  did  not  always  correspond  to  acoustics 
data  collection  intervals  (see  Figure  7).  Thus  a  significant  percentage  of  the  acoustics  data  were  not 
processed  beyond  this  step  due  to  a  lack  of  GPS  positioning  support.  Second,  because  of  a 
downloading  error  aboard  ship,  only  single-frequency  GPS  positioning  from  the  Ashtech  3DF 
receiver  was  available  to  support  calibration  Une  data  reduction  for  day  142.  Third,  since  an  eight- 
point  interpolator  was  used,  some  GPS  position  solutions  on  each  end  of  a  continuous  span  were 
not  available. 

In  addition,  ship  maneuvers  during  any  2()-sec  time  span  and  breaks  in  the  available  GPS 
or  acoustics  data  sets  impacted  the  procedure  used  to  determine  ship  position  at  acoustic  signal 
receive  times.  This  procedure  estimated  the  ship  velocity  vector  from  interpolated  GPS  ship 
positions  at  consecutive  20-sec  pulse  epochs  and  then  extrapolated  ship  position  forward  to  pulse 
receive  times  using  this  estimated  vector.  Later,  this  approach  was  replied  by  the  direct  use  of  the 
Ashtech  3DF  receiver  data. 

The  Ashtech  3DF  receiver  produces,  among  other  quantities,  estimates  of  ship  speed  and 
course-over-ground.  The  latter  quantity  differs  from  estimated  ship  heading  insofar  as  it  relates  to 
the  orientation  of  the  ship  velocity  vector  not  to  the  orientation  of  the  ship  centerline.  These 
quantities  can  differ  by  several  degrees  depending  on  ship  maneuvers,  currents,  and  other 
conditions.  In  the  latter  phases  of  data  processing,  the  3DF  estimates  of  speed  and  course-over- 
ground  were  used  to  extrapolate  ship  position  from  pulse  emit  to  pulse  receive  times  for  each 
measurement  event 


CONNECTING  THE  GPS  ANTENNA  AND  THE  ACOUSTIC  TRANSDUCER 

To  provide  additional  completeness  in  the  observation  equation,  it  was  necessary  to  account 
for  the  fact  that  the  GPS  antenna  to  which  the  ship  position  solutions  refer  was  mounted  several 
decks  above  the  waterline;  whereas,  the  acoustic  transducer  was  mounted  below  the  hull.  Since 
the  ship  was  constantly  undergoing  changes  in  roll,  pitch,  and  heading,  it  was  necessary  to 
determine  the  instantaneous  orientation  of  the  displacement  vector  connecting  these  points  in  the 
WGS  84  reference  frame.  The  GPS  ship  positions  were  then  transferred  to  the  transducer  via 
vector  addition. 

For  this  deployment,  the  ship  dual-frequency  GPS  antenna  was  antenna  1  of  the  four- 
antenna  array  connected  to  the  Ashtech  3DF  attitude  receiver.  Antenna  1  was  located  16  ft  to  the 
port  of  the  ship’s  centerline,  and  5  ft  aft  of  the  transducer  point.  The  transducer  was  located  2  ft 
starboard  of  the  centerline  of  the  ship.  When  deployed,  it  is  extended  4  m  below  the  hull. 
Figures  8  and  9  provide  the  geometry  relating  the  transducer  location  (retracted)  and  GPS  anteimas 
1  and  2.  In  the  ship-fixed  reference  frame,  the  vector  from  antenna  1  to  the  transducer  (retracted 
position)  is  given  by  the  components  (-18.0,  -5.0, 32.5)  in  units  of  feet.  When  the  transducer  was 
deployed,  the  z-component  of  this  vector  was  modified  to  reflect  its  extension  away  from  the  hull. 
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The  Ashtech  3DF  receiver  collects  C/A  single-frequency  GPS  range  and  phase 
observations  using  a  four-antenna  array  and  produces  a  three-dimensional  attitude  history  for  the 
plane  containing  the  antenna  array.  In  this  particular  configuration,  the  antennas  were  deployed  on 
different  decks  of  the  ship  and  were  not  aligned  to  the  ship  centerline.  Therefore,  adjustments 
needed  to  be  applied  to  the  3DF  attitude  results  to  determine  the  roll,  pitch,  and  heading  of  the  ship. 
After  considering  the  relative  locations  of  the  antennas  on  the  ship’s  decks,  it  was  determined  that 
the  pitch,  roll,  and  heading  offsets  of  the  antenna  array  relative  to  the  body-fixed  axes  of  the  ship 
were  -18.6,  -0.5,  and  -3.1  deg,  respectively.  The  3DF  attitude  data  were  corrected  to  reflect  the 
attitude  dynamics  of  the  ship  reference  frame. 

To  determine  the  position  of  the  transducer  in  the  WGS  84  reference  frame,  a  set  of  rotation 
matrices  had  to  be  applied  to  transform  the  ship-fixed  displacement  vector  between  antenna  1  and 
the  transducer  into  the  WGS  84  frame.  This  rotation  was  accomplished  using  the  following 
transformation  equation. 


XwGS  84  =  Rz(-(  1 80+A,))Ry(90-<p)Rz(90-Hj))Rx(-(pp)Ry(-\jr)Xship  ( 1 ) 

where 

XwGS  84  =  (x,y,z)wGS  84  (2) 

and  Xship  me  the  Cartesian  coordinates  of  the  antenna-transducer  displacement  vector  in  the  body- 
fixed  (ship)  frame.  This  transformation  equation  consists  of  a  rotation  into  the  north-east-up  frame 
followed  by  a  rotation  into  the  geodetic  reference  frame.  The  rotation  angles  \|r,  <pp,  <)>,  (p,  and  X 
used  in  this  transformation  are  the  roll,  pitch,  heading,  latitude,  and  longitude,  respectively,  of  the 
ship.  All  angular  quantities  are  derived  from  the  GPS  observables.  The  rotation  matrices  Rx,  Ry, 

Rz  are  the  standard  Euler  angle  rotations.^ 

Once  the  displacement  vector  is  known  in  WGS  84,  it  is  added  to  the  WGS  84  Cartesian 
coordinates  of  antenna  1,  determined  using  dual-frequency  GPS  observables,  to  yield  the  WGS  84 
Cartesian  coordinates  of  the  transducer  at  emit  or  receive  times. 


EDITING  TO  PRODUCE  COMMON  DATA  INTERVALS 

Figure  7  shows  that  the  data  collection  interval  during  this  experiment  was  not  consistent 
across  data  types.  As  various  phases  of  acoustics  data  processing  were  performed,  as  discussed  in 
detail,  the  volume  of  the  acoustics  data  considered  varied  depending  on  the  degree  to  which  it  could 
be  supported  by  the  availablity  of  other  data  types.  All  data  sets  were  edited  at  each  step  in  the 
processing  to  produce  the  largest  common  data  interval  for  that  processing  phase. 
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,  INITIAL  ADJUSTMENT  RESULTS 


A  first  adjustment  using  only  acoustics  data  was  performed  incorporating  all  data  files  with 
the  exception  of  hour  0  on  day  143.  This  one  data  set  was  not  initially  processed  because  of  some 
difficulty  with  reading  the  file  on  the  CDC  995  system.  A  total  of  over  7000  acoustic  observations 
were  processed.  These  data  were  acquired  from  the  three  USGS  transponders  (T5,  T6,  T7)  and 
one  NOAA  transponder  T4  (see  Figure  5).  Initial  coordinates  for  transponder  and  transducer 
positions  were  provided  by  the  USGS  in  a  local  reference  frame  in  which  the  NOAA  transponder 
coordinates  were  known.  The  adjustment  assigned  an  initial  uncertainty  of  5  m  to  each  acoustic 
observation  and  a  position  uncertainty  of  20  m  to  each  ship  position  coordinate.  The  ship  tracks 
during  this  calibration  line  collection  period  are  shown  in  Figure  10. 
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RGURE  10.  CALIBRATION  LINE  SHIP  TRACKS  OVER  USGS  TRANSPONDER 
ARRAY  ON  DAY  142,  HOUR  18  TO  DAY  143,  HOUR  3 


The  purpose  of  this  initial  adjustment  was  to  examine  the  internal  consistency  of  the  data 
set,  to  establish  a  basis  for  data  editing  using  residuals  of  fit,  and  to  determine  appropriate  data 
weighting  for  subsequent  adjustments.  The  starting  transponder  coordinates  and  the  resulting 
adjustments  to  these  coordinates  are  given  in  Table  2. 


I 
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TABLE  2.  STARTING  TRANSPONDER  COORDINATES*  AND  RESULTS  FROM  INITIAL  ADJUSTMENT 

(units  in  meters) 


Unit 

X 

Y 

Z 

AX 

AY 

AZ 

T5 

11974 

14137 

2205 

-2.0 

-4.6 

7.3 

T6 

12902 

14627 

2205 

-1.8 

5.7 

11.3 

T7 

12240 

14800 

2205 

8.3 

-2.4 

10.3 

T4 

15443 

18396 

1987 

8.8 

-4.7 

31.8 

*LocaI  reference  frame 


Several  assumptions  were  made  at  this  stage  in  the  processing  and  important  corrections  to 
the  data  had  not  been  applied.  For  instance,  acoustic  refraction  corrections  could  not  yet  be 
computed  because  the  necessary  sound  velocity  profile  was  not  available  from  NFS,  therefore  a 
constant  sound  velocity  of  1482.62  m/sec  was  used  for  data  processing  initially.  Also,  a 
transponder  internal  electronic  delay  of  0.015  sec  was  adopted  for  each  instrument  as  suggested  by 
USGS.  These  factors  are  the  most  likely  causes  for  the  z-coordinate  (downward)  adjustment  in 
Table  2  for  all  transponders  units. 

Figure  1 1  illustrates  the  adjustments  occurring  to  the  initial  coordinates  for  the  transponders 
as  the  normal  equations  were  sequentially  merged  and  solved.  This  scheme  was  followed 
beginning  with  the  normal  equations  for  calibration  track  18  (hour  18,  day  142)  through  track  03 
(hour  3,  day  143).  Figure  1 1  shows  the  results  at  several  stages  of  this  process;  the  last  step 
18-03  being  the  solution  based  on  all  acoustics  data  with  the  exception  noted  above.  For  the 
coordinates  of  transponders  T5,  T6,  and  T7  these  sequential  adjustments  were  more  consistent 
(clustered)  than  for  the  mote  distant  NOAA  unit  T4. 

Figure  12,  corresponding  to  Figure  11,  provides  the  coordinate  covariance  evolution  as 
track  data  were  merged.  Although  not  necessarily  a  valid  representation  of  the  statistical 
uncertainty  for  the  solution,  these  numerical  values  do  provide  some  general  insight  into  the  data 
set’s  ability  to  provide  accurate  positioning:  first,  the  data  set  supports  a  stronger  solution  for  depth 
than  it  does  for  the  horizontal  coordinates  of  the  USGS  array;  and  second,  the  ability  to  position 
the  more  distant  NOAA  transponder  T4  is  much  weaker  due  to  the  actual  geometry  of  the  data 
acquisition. 
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HGURE  11.  INITIAL  ACOUSTIC  DATA  SOLUTION-COORDINATE 
ADJUSTMENTS  AS  CALIBRATION  LINE  DATA  ARE  MERGED 


TRACKS 


TRACKS  18-21 
TRACKS  18-23 

TRACKS  18-03 


COORDINATE 


HGURE  12.  COORDINATE  COVARIANCE  AS  TRACK  DATA  ARE  MERGED 
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This  adjustment  also  provided  a  significant  view  into  the  nature  of  the  acoustics  data. 
Figures  13  and  14,  for  example,  are  plots  of  the  post-adjustment  data  residuals  for  two  of  the 
calibration  tracks.  Although  there  are  systematic  trends  in  the  residuals  evident  for  each 
transponder,  it  can  readily  be  seen  that  spurious  data  points  and  large  systematic  errors  are  present. 
Each  of  the  acoustic  files  was  then  edited  to  eliminate  suspect  data.  Figures  15  and  16  are  two 
examples  of  the  residuals  of  fit  after  data  editing,  prior  to  a  second  adjustment  using  edited  acoustic 
data  sets. 


SECONDS  FROM  START  OF  HOUR 

HGURE  13.  POST-ADJUSTMENT  RESIDUALS  (HOUR  19) 


SECONDS  FROM  START  OF  HOUR 


FIGURE  14.  POST-ADJUSTMENT  RESIDUALS  (HOUR  23) 
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HGURE  15.  EDITED  POST-ADJUSTMENT  RESIDUALS  (HOUR  19) 
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HGURE  17.  STANDARD  DEVIATIONS  OF  HOURLY  ACOUSTICS 
DATA  SETS  BEFORE  AND  AFTER  EDITING 


The  mean  and  standard  deviation  of  each  hourly  data  set  were  computed  prior  to  and  after 
data  editing  using  values  for  the  coordinates  from  the  initial  adjustment  Table  3  provides  these 
results.  The  data  set  multiplier  given  in  the  Table  3  was  the  factor  by  which  the  standard  deviation 
of  the  pre-edited  data  was  multiplied  to  establish  a  data  rejection  (editing)  criteria.  In  general,  it 
appears  that  the  low-frequency  acoustics  data,  after  being  converted  to  two-way  ranges  using  a 
fixed  sound  velocity,  have  standard  errors  of  from  5  to  12  m  prior  to  correcting  for  other  factors 
such  as  anomalous  refraction  effects. 


TABLE  3.  HOURLY  DATA  SET  STATISTICAL  MEASURES,  PRE-  AND  POST-EDITED 

(units  in  meters) 


Hour 

Mean 

Sigma 

Edited 

Mean 

Edited 

Siema 

Multiplier 

18 

3.9 

34.3 

0.1 

11.5 

3.0 

19 

-1.1 

12.4 

-0.8 

8.5 

2.5 

20 

2.0 

11.0 

0.7 

4.9 

2.5 

21. 

-1.3 

8.7 

-1.4 

6.1 

2.5 

22 

-1.3 

10.7 

-0.8 

5.6 

2.5 

23 

0.0 

17.5 

-1.9 

6.3 

1.0 

01 

-1.5 

12.0 

-1.9 

6.7 

1.5 

02 

-0.1 

9.4 

-0.8 

7.2 

2.5 

03 

-0.2 

16.8 

-0.5 

8.0 

1.5 

18 


NSWCDD/TR-93/643 


A  second  adjustment  was  then  made  using  these  edited  files.  This  adjustment  was 
performed  using  exactly  the  same  conditions  as  the  first  but  included  only  the  edited  data  sets. 
Table  4  shows  the  impact  of  data  editing  on  the  parameter  adjustments.  Significant  changes  to 
some  of  the  parameters  resulted,  implying  that  data  editing  would  be  critical  to  obtaining  a  final 
solution. 


TABLE  4.  RESULTS  FROM  INITIAL  ADJUSTMENT  AND 
ADJUSTMENT  USING  EDITED  ACOUSTICS  DATA 
(units  in  meters) 


Initial  Adjustment 

Adjustment  Using  Edited  Data 

Unit 

AX 

AY 

AZ 

AX 

AY 

AZ 

T5 

-2.0 

-4.6 

7.3 

-3.1 

4.5 

1.4 

T6 

-1.8 

5.7 

11.3 

-5.1 

4.9 

3.5 

T7 

8.3 

-2.4 

10.3 

10.2 

-9.2 

3.9 

T4 

8.8 

-4.7 

31.8 

-13.3 

-4.4 

26.7 

ADJUSTMENT  INTO  THE  WGS  84  REFERENCE  FRAME 


Even  though  these  initial  solutions  provided  good  insight  into  the  characteristics  of  the 
acoustics  data,  they  are  of  little  value  to  those  interested  in  geodynamics  or  geodetic  positioning 
because  they  are  not  linked  to  any  global  reference  frame.  Ideally,  the  positions  of  these 
instruments  would  be  known  with  high  accuracy  in  a  well-established  global  reference  frame  so 
that,  in  the  future,  any  detected  spreading  rates  derived  from  high-frequency  acoustics  could  be 
expressible  as  a  vector  quantity  in  the  same  frame.  This  would  allow  a  more  comprehensive 
merging  of  such  marine  experimental  results  with  terrestrial  geodynamical  findings.  This  is  where 
GPS  plays  a  critical  role. 

As  described,  the  USGS  Juan  de  Fuca  ridge  experiment  was  supported  by  two  dual¬ 
frequency  GPS  receivers,  one  antenna  located  on  an  upper  deck  of  the  iMney  Chouest  and  a 
second  located  in  Newport,  Oregon  at  a  point  of  known  geodetic  control  in  WGS  84.  The  data 
collection  at  the  two  receivers  allowed  the  ship’s  track  to  be  determined  in  the  WGS  84  reference 
>  frame  with  an  accuracy  sufficient  for  processing  the  acoustics  data  set.  Since  the  ship  positions 

I  derived  from  GPS  are  in  the  WGS  84  reference  frame,  their  use  in  the  acoustics  data  processing  as 

weighted  parameters  ensures  that  the  transponder  coordinates  are  determined  within  the  same 
geodetic  system. 

At  this  stage  in  the  data  analysis,  two  additional  solutions  were  made  supported  by  GPS- 
I  derived  ship  positions.  The  first  (I)  included  all  acoustics  data  acquired  on  day  142  beginning  at 

hour  18  and  ending  on  day  143  at  hour  3  with  the  exception  of  hour  0.  The  second  (II)  included 
only  the  data  used  in  the  first  solution  acquired  during  day  143,  hour  1  to  hour  3.  Again,  as  in  the 

[acoustics  data  solutions  described  earlier,  refraction  corrections  were  not  applied.  In  addition,  the 
relationship  between  the  ship’s  GPS  antennas  and  the  transducer’s  position  was  purposefully 
neglected  at  this  stage,  since  the  Ashtech  3DF  data  for  day  142  had  not  yet  been  received.  Both 
solutions  used  the  edited  acoustics  data  (edited  to  remove  spurious  observations  and  acoustics  data 
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when  GPS  solutions  were  not  available  (see  Figure  7)).  Each  solution  used  data  weighting 
consistent  with  sample  statistics  developed  from  examining  previous  residuals  of  fit. 

These  solutions  were  processed  assuming  that  the  GPS  ship  positions  were  accurate  to 
20  m.  As  it  turned  out,  on  the  ship  at  the  end  of  day  142,  the  GPS  data  was  downloaded  to  tape. 
However,  instead  of  downloading  dual-frequency  data,  only  single-frequency  Ashtech  3DF  data 
were  recorded  on  tape.  Both  of  the  GPS  ship-based  receivers  were  connected  to  antenna  1  of  the 
Ashtech  3DF  system.  This  mistake  resulted  in  the  loss  of  dual-frequency  data  on  day  142  so  that 
the  first  GPS/acoustics  data  solution  had  a  mix  of  GPS  single-frequency  and  dual-frequency 
determined  ship  positions.  After  processing  the  acoustics  data  on  a  hourly  basis,  the  normd 
equations  were  appropriately  merged  to  produce  the  two  solutions  in  the  WGS  84  reference  frame. 
The  second  solution  based  on  dual-frequency  GPS  contained  significantly  less  data. 

To  compare  these  GPS-based  solutions  with  the  acoustics  (edited  data)  solution  in  a  local 
reference  frame,  baseline  distances  were  computed  in  each  frame.  The  invariant  nature  of  this 
quantity  with  respect  to  reference  frame  makes  it  useful  for  this  purpose.  Figure  18  provides  the 
baseline  lengths  from  these  three  solutions.  Figure  19  provides  the  differences  in  these  baseline 
lengths.  Table  5  summarizes  these  and  other  solutions  performed  during  this  study. 


BASEUNE 

HGURE  18.  BASELINE  LENGTHS  FROM  PRELIMINARY  SOLUTIONS 


A 

I 

II 


20 


NSWCDD/TR-93/643 


80 


T5-T6  T5-T7  T5-T4  T6-T7  T6-T4  T7-T4 


BASELINE 

HGURE  19.  BASELINE  LENGTH  DIFFERENCES  FROM  PRELIMINARY  GPS 
SOLUTIONS  WHEN  COMPARED  TO  ACOUSTICS-ONLY  SOLUTION  (A) 


TABLES.  SOLUTION  SET DEHNinONS:  TRANSPONDERS  T5, T6,  T7,  T4 
(DAY  142,  HOUR  18  TO  DAY  143,  HOUR  3) 


Designator 

Acoustics 

GPS  Position 

GPS 

Attinide 

Refraction 

CorrectkM 

Transmit  Time 

Receive  Time 

A 

18-03 

No  Hour  (X) 

No 

No 

No 

No 

I 

18-03 

No  Hour  00 

Day  142  LI  Only 
Dav  143  LI  &  L2 

Interpolate  using  GPS 
20-sec  solutions 

No 

No 

n 

01-03 

LI  &  L2  Solution 

Interpolate  using  GPS 
20-sec  solutions 

No 

No 

m 

01-03 

LI  &  L2  Solution 

Interpolate  using  GPS 
20-S6C  solutions 

Yes 

No 

IV 

01-03 

LI  &  L2  Solution 

Use  GPS-derived 
speed  and  course 

Yes 

No 

V 

00-03 

LI  &  L2  Solution 

Use  GPS-derived 
speed  and  course 

Yes 

No 

VIA 

LI  &  L2  Solution 

Use  GPS-derived 
speed  and  course 

Yes 

Yes 

VI B* 

00-03 

LI  &  L2  Solution 

Use  GPS-derived 
speed  and  course 

Yes 

Yes 

"'Solutioii  includes  T5,  T6,  and  T7  coordinates  only. 
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The  results  provided  in  Figure  19  imply  that  the  use  of  single-frequency  GPS  differential 
positions  for  the  ship  at  a  range  of  300  nmi  may  not  prove  to  be  satisfactory.  Although  the  amount 
of  acoustics  data  present  in  the  first  GPS  solution  is  greater  than  the  second  solution,  the  average 
baseline  difference  was  45.4  m  as  compared  to  15.2  m  resulting  from  less  data  but  supported  by 
dual-frequency  GPS  differential  positioning.  Notice  that  in  each,  the  scale  of  the  solution  is 
excessive  when  compared  to  the  acoustics  solution  in  the  local  reference  frame. 


SHIP  ATTITUDE  DATA 


In  the  previous  solutions,  the  relationship  between  the  GPS  antenna  1  cormected  to  both  the 
Ashtech  3DF  and  LDXII  receivers  and  the  hull  mounted  transducer  had  been  ignored.  During  the 
next  stage  of  data  processing,  the  3DF  attitude  results  providing  the  real-time  geometry  of  the 
antenna-to-transducer  displacement  vector  were  incorporated  into  the  observation  equation. 

The  edited  acoustic  data  files  on  day  143,  hours  1-3,  containing  dual-frequency  GPS  ship 
positions,  were  augmented  with  3DF-derived  values  of  ship  heading,  pitch,  and  roll  at  each 
transducer  emit  time.  The  vector  between  antenna  1  and  the  transducer,  determined  in  the  ship 
body-fixed  system,  was  transformed  using  equation  (1)  into  the  WGS  84  reference  frame.  This 
displacement  vector  was  then  added  to  the  position  vector  of  antenna  1  derived  from  GPS 
measurements  to  determine  the  WGS  84  position  vector  components  for  the  transducer  at  the  pulse 
emit  time. 

As  in  previous  solutions,  ship  positions  at  the  receive  times  of  the  returning  pulses  were 
interpolated  from  GPS-derived  ship  positions  spaced  at  20  sec.  In  addition,  the  ship  attitude  was 
considered  constant  over  the  3-  to  6-sec  time  intervals  following  pulse  emit  This  simplifying 
assumption  is  generally  supported  by  the  data. 

After  the  data  were  processed,  the  distances  between  transponder  pairs  were  recomputed 
using  the  estimated  WGS  84  coordinates  and  then  differenced  with  the  distances  obtained  from  the 
acoustics-only  solution.  Figure  20  presents  these  distance  differences  as  well  as  repeats  the  results 
presented  in  Figure  19.  Although  Ae  modeling  (3DF  data)  incorporated  into  this  last  solution  (HI) 
was  expected  to  have  enhanced  the  previous  solution  (designated  as  GPS  143),  the  mean 
difference  in  the  baseline  lengths  actually  increased  from  15.2  to  16.8  m. 
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80 


TS-T6  T5-T7  T5-T4  T6-T7  T6-T4  T7-T4 

BASELINE 

HGURE  20.  BASELINE  LENGTH  DIFFERENCES  FOR  SEVERAL  PRELIMINARY  SOLUTIONS 


As  a  further  step  toward  establishing  a  convergence  between  the  acoustic  data-only  results 
and  those  supported  by  GPS,  the  procedure  used  to  determine  GPS  ship  coordinates  at  times  of 
acoustic  signal  return  was  reexamined.  Although  several  options  for  changing  the  approach 
existed,  it  was  decided  that  the  ship’s  speed  and  course-over-ground  from  the  Ashtech  3DF 
receiver  would  be  used.  This  approach  would  still  minimize  the  number  of  files  being  processed, 
but  would  circumvent  some  of  the  weaknesses  inherent  with  interpolating  GPS  20-sec  positions. 

Since  speed  and  course-over-ground  were  referenced  to  a  local  level  frame,  it  was 
necessary  to  transform  these  quantities  into  the  WGS  84  reference  frame  as  was  previously  done 
for  the  3DF  attitude  data.  Then,  assuming  that  these  quantities  were  constant  over  the  next  3  to 
6  sec,  the  ship’s  displacement  vector  over  each  pulse  emit/retum  interval  was  calculated. 
Transducer  coordinates  at  receive  times  were  then  determined  in  WGS  84  and  used  in  the  acoustic 
data  processing. 

Figure  21  adds  to  the  previous  baseline  comparisons  the  distance  differences  between  this 
next  solution  (IV)  incorporating  3DF  attitude,  speed,  and  course-over-ground  and  the  acoustics- 
only  solution.  The  mean  baseline  difference  was  now  6.9  m;  however,  the  T6-T7  difference 
remained  large. 
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HGURE  21.  BASELINE  LENGTH  DIFFERENCES  FROM  ACOUSTICS-DATA-ONLY  SOLUTION 


Next,  the  observations  from  day  143,  hour  0  were  added  to  the  data  processing.  These 
acoustics  observations  were  supported  by  dual-frequency  GPS  positioning  but  were  previously  not 
processed  because  of  difficulty  in  reading  the  file.  The  results  of  the  baseline  length  comparisons 
from  this  solution  (IV)  with  the  acoustics-data-only  solution  are  also  given  in  Figure  21.  The  mean 
difference  in  baseline  length  was  2.9  m,  with  the  T6-T7  baseline  difference  still  excessive.  The 
ship  track  pattern  showing  where  this  data  (day  143,  hour  0  to  3)  was  acquired  is  given  in 
Figure  22,  which  can  be  contrasted  with  Figure  10. 
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HGURE  22.  CALIBRATION  TRACKS  FOR  DAY  143,  HOUR  0-3 


ACOUSTIC  REFRACTION  MODELING 


As  a  final  processing  step,  refraction  corrections  were  estimated  and  applied  to  the 
acoustics  data  using  a  layered  model  for  the  sound  velocity  profile  within  the  water  column.  This 
profile  was  developed  from  expendible  bathythermograph  observations  collected  in  1986  by 
NOAA  in  the  vicinity  of  the  USGS  experiment  The  profile  is  given  in  Table  6. 

The  acoustic  pulse  reaching  each  transponder  arrives  along  a  refracted  path  that  depends  on 
the  temperature,  pressure,  and  salinity  within  Ae  water  column.  These  variations  are  observed  as  a 
function  of  depth  and  are  converted  into  a  layered  model  parameterized  by  depth  and  sound 
velocity.  Such  layered  profiles  (see  Figure  23)  support  ray  tracing  based  on  Shell’s  Law.^ 
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TABLE  6.  1986  NOAA  SOUND  VELCX:iTY  PROFILE 
IN  WATER  COLUMN  NEAR  JUAN  DE  FUCA  RIDGE 
EXPERIMENT  AREA 


■KlSHi 

Velocity 

(ft/sec) 

0.0 

4954.40 

3.281 

4954.40 

32.808 

4953.41 

65.617 

4952.43 

98.425 

4950.46 

131.234 

4936.35 

164.042 

4915.68 

246.063 

4888.78 

328.084 

4882.87 

492.126 

4863.19 

656.168 

4863.85 

984.252 

4853.67 

1312.336 

4846.13 

1640.420 

4843.83 

2460.630 

4848.75 

3280.840 

4855.31 

4101.050 

4862.86 

4921.260 

4871.06 

5741.470 

4879.59 

6561.680 

4890.09 

7381.890 

4903.54 

Velocity 

profile 


nGURE23.  REFRACTION  IN  A  LAYERED  MEDIUM 
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As  previously  discussed,  measurements  within  the  water  column  were  acquired  during  the 
deployment  of  the  USGS  transponders  by  instrumentation  mounted  onboard  the  deep-submersible 
Turtle.  However,  difficulty  with  software  prevented  the  actual  sound  velocity  profile  from  being 
available  in  time  to  support  this  data  analysis.  The  NOAA  profile  given  in  Table  6  was  used  in  the 
interim. 


Each  acoustic  observation  was  corrected  for  refraction  using  the  following  procedure: 
signal  emit  and  receive  times  were  used  to  determine  positions  for  the  transducer  in  the  WGS  84 
reference  frame  based  on  GPS.  Having  approximate  coordinates  for  the  transponders  in  the  same 
frame,  as  determined  through  previous  processing,  ray  tracing  through  the  profile  in  Table  6  was 
used  to  determine  the  amount  of  time  the  acoustic  pulse  took  to  reach  then  return  from  each 
transponder.  Using  approximate  geometry,  the  angle  of  depression  at  the  transducer  to  each 
transponder  was  determined  to  initiate  the  ray  trace.  The  tracing  terminated  at  the  bottom  depth  of 
each  transponder  at  an  offset  position  from  each  instrument.  These  horizontal  offsets  were  used  to 
refine  the  depression  angles  at  the  surface;  then  ray  tracing  was  repeated.  This  process  was  iterated 
until  the  traced  rays  reached  each  transponder  location  to  within  a  meter.  The  process  usually 
converged  within  three  iterations  for  transponders  T5,  T6,  and  T7.  For  the  more  distant 
instrument  T4,  the  process  required  additional  care  in  adjusting  the  depression  angle  at  each 
iteration  and  took  around  10  iterations  to  converge. 

Once  the  ray  paths  to  and  from  the  transponder  were  determined,  the  time  of  transit  along 
the  bent  paths  was  calculated  and  differenced  with  an  estimated  time  of  transit  along  the  straight- 
line  geometric  paths  that  would  have  been  followed  had  the  sound  velocity  profile  been  constant. 
The  constant  velocity  used  in  this  calculation  was  a  weighted  average  determined  for  each 
transponder  using  the  NOAA  profile.  The  weights  were  based  on  the  thickness  of  each  profile 
layer  with  the  bottom  depth  of  each  instrument  terminating  the  profile  in  each  calculation.  The 
difference  between  the  round-trip  travel  time  along  the  bent  rays  and  along  the  straight-line  paths 
was  used  as  an  estimate  of  the  refraction  delay  for  each  data  point  and  was  applied  as  a  correction. 

This  approach  for  correcting  for  the  effects  of  water  column  refraction  was  certainly  not 
optimal  since  it  incorporated  certain  assumptions  and  a  NOAA  profile  that  was  only  approximate. 
An  alternate  procedure  for  modeling  the  observation  equation  might  circumvent  some  weaknesses 
in  the  assumptions.  This  alternate  approach  is  summarized  in  the  appendix. 

Although  not  optimal,  these  refraction  corrections  were  applied  to  the  acoustics  data 
available  during  day  143  when  dual-frequency  GPS  observations  supported  ship  positioning. 
Figures  24  and  25  provide  one  example  of  the  effect  that  this  refraction  correction  scheme  had  on 
the  pre-adjustment  data  residuals.  Figure  24  is  a  graph  of  residuals  for  transponder  T5  on  day 
143,  hour  1,  prior  to  correction  for  refraction.  Figure  25  provides  the  data  residuals  after 
refraction  correction.  It  is  evident  that  the  refraction  correction  scheme  has  removed  much  of  the 
high  frequency  noise  on  the  acoustic  data  due  to  refraction  effects.  However,  systematic  effects 
are  still  present  in  the  data.  For  instance,  the  step  in  the  residuals  occurring  near  data  point  55  can 
be  associated  with  a  change  in  the  GPS  satellites  being  tracked.  This  step  was  also  present  in  the 
residuals  from  hour  1  for  the  other  transponders.  It  is  suspected  that  much  of  the  remaining 
systematic  trends  in  the  residuals  are  due  to  GPS  positioning  errors  and  inadequacies  in  the 
refraction  model. 
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HGURE  24.  TRANSPONDER  T5  PRE-ADJUSTMENT  RESIDUALS 
WITH  NO  REFRACTION  CORRECTION  (DAY  143,  HOUR  1) 


HGURE  25.  TRANSPONDER  T5  PRE-ADJUSTMENT  RESIDUALS 
WITH  REFRACnON  CORRECTION  (DAY  143,  HOUR  1) 


After  incorporating  refraction  correction  into  the  adjustment  model,  data  for  day  143  was 
processed  to  determine  the  coordinates  of  the  transponder  network  in  the  WGS  84  reference  frame. 
All  data  having  pre-adjustment  residuals  greater  than  50  m  were  edited  from  the  solution.  An  initial 
solution  using  this  data  included  coordinate  corrections  for  all  four  transponders.  After  data 
processing,  residuals  of  Ht  and  all  combinations  of  baseline  lengths  were  computed.  Figure  26 
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illustrates  the  residuals  after  this  adjustment  for  hour  1.  Two  points  need  to  be  made  regarding 
Figure  26.  First,  systematic  trends  in  the  residuals  remain,  and  second,  numerous  data  residuals 
for  the  NOAA  transponder  T4  are  aggregated  (see  box  in  figure)  away  from  those  for  transponders 
T5,  T6,  and  T7. 


DATA  POINTS 


HGURE  26.  POST-ADJUSTMENT  RESIDUALS  FOR  ALL 
TRANSPONDERS  (DAY  143,  HOUR  1) 


Considering  the  data  acquisition  geometry  presented  in  Figure  22  and  the  residuals 
associated  with  transponder  T4,  a  second  solution  was  made  deleting  all  T4  data  from  the 
processing.  It  was  known  that  the  geometry  (figure  of  merit)  for  the  solution  for  the  coordinates 
for  T4  was  much  weaker  and,  because  of  the  path  lengths,  the  refraction  corrections  may  be 
somewhat  (more)  suspect. 

These  last  two  solutions  were  designated  VIA  and  VIB,  respectively.  Figure  27  depicts  the 
resulting  baseline  length  differences  between  several  of  the  solutions  discussed  to  this  point  using 
data  supported  by  dual-frequency  GPS  tracking,  GPS  attitude,  speed,  and  course-over-ground 
(see  Table  5  for  solution  definition  summary). 
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HGURE  27.  BASELINE  DIFFERENCES  BETWEEN  VARIOUS  GPS- 
SUPPORTED  SOLUTIONS  AND  ACOUSTICS-ONLY  SOLUTION 


A  comparison  of  the  baseline  solutions  from  V  with  VIA  shows  that  the  inclusion  of 
refraction  correction  has  an  impact  on  the  resulting  lengths  ranging  from  less  than  1  m  ((T5-T6) 
and  (T7-T4))  to  over  15  m  (T6-T7).  When  a  second  solution  with  refraction  correction  was 
performed  with  T4  data  removed,  the  change  in  the  baselines  between  T5,  T6,  and  T7  ranged  from 
2  to  10  m. 


FINAL  RESULTS  AND  CONCLUSIONS 


To  date,  a  number  of  solutions  have  been  attempted,  which  incorporated  differences  in 
modeling  and  variations  in  the  data  sets.  The  solutions  that  combined  GPS  observables  and 
acoustic  measurements  have  been  compared  to  an  acoustics-data-only  solution  that  determined  ship 
and  transponder  positions  in  a  local  reference  frame.  Comparisons  of  different  solutions  were 
made  using  distances  between  transponders  computed  from  estimated  coordinates.  Since  these 
quantities  are  invariant  with  respect  to  coordinate  system,  they  offer  one  method  of  determining 
consistency  between  solutions.  Figure  28  provides  the  differences  in  the  distances  between 
transponders  based  on  comparing  the  acoustics-data-only  solution  with  seven  other  solutions 
derived  using  GPS  and  acoustics  data.  These  latter  solutions  are  in  WGS  84  and  were  defined  in 
Table  5.  The  six  baseline  combinations  shown  in  Figure  28  were  computed  for  all  GPS-based 


30 


NSWCDD/TR-93/643 


solutions  with  the  exception  of  VIB,  which  did  not  include  data  from  transponder  T4.  All  data 
reductions  using  GPS  were  based  on  a  constant  sound  velocity  in  the  water  column.  For  solutions 
VIA  and  VIB  that  constant  was  computed  for  each  transponder  using  a  1986  NOAA  sound  velocity 
profile. 
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T5-T6  T5-T7  T5-T4  T6-T7  T6-T4  T7-T4 

BASELINE 


HGURE  28.  DIFFERENCES  IN  BASELINE  LENGTHS  WHEN  VARIOUS  GPS 
SOLUTIONS  ARE  COMPARED  TO  ACOUSTICS-ONLY  SOLUTION 


Solution  V  incorporated  the  full  use  of  GPS  data  and,  on  average,  agreed  to  2.9  m  with 
baseline  lengths  from  the  acoustics-data-only  solution.  However,  due  to  the  non-availability  of 
dual  frequency  GPS  data  on  day  142,  the  number  of  observations  in  solution  V  was  significantly 
less  than  for  the  acoustics  solution,  which  might  account  for  the  larger  difference  in  the  T6-T7 
baseline  length. 

With  the  addition  of  refraction  correction,  it  was  shown  that  the  baseline  scale  increased  for 
the  USGS  network  but  decreased  for  baselines  extending  to  T4.  As  this  seemed  suspect,  a  final 
solution  VIB  for  T5,  T6,  and  T7  coordinates  was  computed  and  is  given  in  Table  7.  TTie  accuracy 
of  this  solution  remains  to  be  verified  using  directly  measured  acoustic  distances  that  should  be 
retrieved  by  the  USGS  on  a  subsequent  site  visit 
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TABLE  7.  USGS  TRANSPONDER  COORDINATES  IN  WGS  84  FROM  SOLUTION  VIB 

(units  in  meters) 


Unit 

X 

Y 

Z 

T5 

-2941651.8 

-3460767.0 

4459469.7 

T6 

-2940709.0 

-3461114.1 

4459811.7 

T7 

-2941146.4 

-3460589.4 

4459938.1 

RECOMMENDATIONS 


The  full  potential  for  developing  highly  accurate  geodetic  coordinates  from  the  USGS 
experiment  was  not  realized  for  sever^  reasons.  Some  of  these  were  related  to  the  data  acquisition 
phase  of  the  experimenL  while  others  were  a  result  of  limitations  imposed  during  data  processing. 

With  respect  to  data  acquisition,  the  following  limitations  were  present  in  the  data  set. 

•  Simultaneous  collection  of  acoustic  and  GPS  support  data  did  not  always  occur, 
resulting  in  the  non-processing  of  significant  amounts  of  acoustics  data. 

•  GPS  dual-frequency  data  on  day  142  was  lost. 

•  High-frequency  acoustics  measurements  were  commanded  but  not  retrieved, 
preventing  validation  of  geodetic  positioning  results. 

•  Data  acquisition  geometry  (ship  track  pattern)  during  the  calibration  period  did 
not  yield  a  best  figure  of  merit. 

•  An  in  situ  water  velocity  profile  was  never  provided  to  support  data  analysis. 

Data  analysis  was  limited  by  several  of  these  factors  and  by  the  following: 

•  The  processing  technique  to  obtain  GPS  ship  positions  was  not  optimal  and 
limited  the  accuracy  of  the  results. 

•  Acoustics  data  were  processed  using  an  explicit  representation  for  the 
observation  equation  that  mapped  transit  time  into  distance  using  an  assumed 
constant  for  the  sound  velocity. 

•  Acoustics  data  were  corrected  for  refraction  using  a  1986  NOAA  model. 

Based  on  the  circumstances  surrounding  data  collection  and  processing  limitations,  it  is 
recommended  that  the  original  data  set  be  reprocessed  at  a  later  date  using  the  actual  water  velocity 
profile  and  a  more  accurate  GPS  ship  positioning  scheme.  In  addition,  it  is  recommended  that 
during  the  next  site  visit  that  a  similar  data  set  be  acquired  in  a  manner  that  circumvents  some  of  the 
previous  data  acquisition  limitations. 
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ACOUSTIC  OBSERVATION  EQUATIONS 


To  first  order  the  round-trip  travel  time  T  of  the  acoustic  signal  from  the  ship  hull 
transducer  to  transponder  Ti  is  modelled  by  the  explicit  equation 


Ti  =  [R(S.Ti)  +  R(Ti,Ri)]/Vw  +  A  Ti  +  A  Tr  i  =  1 ,2,3,4  (A- 1) 

where  is  the  average  sound  velocity  in  seawater;  R(S.Ti)  is  the  geometric  slant  range  from  the 
transducer  to  transponder  Ti  at  the  emit  time;  R(Ti.Ri)  is  the  geometric  slant  range  from  the 
transponder  to  the  transducer  at  the  time  of  signal  receipt;  ATi  is  the  calibrated  transponder 
electronic  delay:  and  A  7^  is  a  conection  due  to  refraction  in  the  water  column. 


ESTIMATION  EQUATIONS 


Acoustics  data  are  processed  sequentially,  on  a  track-line  basis.  Transducer  positions  at 
pulse  transmit  and  receive  times  are  assumed  to  be  known  with  an  a  priori  uncertainty.  The 

standard  error  on  these  ship  positions  Os  or  orj  varies  with  the  methods  used  to  position  the  ship, 
geometry,  and  other  factors. 

For  a  given  track  of  data,  observations  are  grouped  as  discrete  events.  Each  event  consists 
of  the  data  associated  with  a  unique  transmitted  acoustic  pulse  and  with  the  positions  of  up  to  four 
transponders  and  five  ship  locations.  Each  event  is  processed  separately. 

Given  an  event,  the  observation  equation  (A-1)  is  linearized  about  the  a  priori  coordinates 
for  the  transponder  array  and  about  coordinates  for  the  transducer  positions  (emit  and  receive): 

V  =  AX  +  L  (A-2) 

Least  squares  normal  equations  are  formed  for  the  event,  which  include  coordinate  corrections  for 
transponders  Ti  (i  =  1, 2,3,4)  and  ship  positions  S  and  Ri  (i  =  1, 2,3,4).  Separating  transponder 
and  ship  position  parameters,  the  normal  equations  are  given  in  partitioned  form  as 

(AtWA)X  +  AtWL  =  0 


NX  +  U=0 


I 

I 

I 


or 
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where 


A  =  (Ai,  A2] 


L  —  Lo  ■  L-b 

Xt  =  [XT,  Xs  Xr]  =  [Xi,X2] 


N  = 


Nil  Ni2 

LN21  N22J 


U  = 


Ui 

LU2, 


Nii=Ai«WAi 
Ni2  =  A1WA2 
N21  =  Ni2^ 
N22  =  A2*WAi 


(A-3) 


Ui  =  Ai^WL 
U2  =  A2‘WL 

The  design  matrix  A  is  developed  by  differentiating  the  observation  equations  with  respect  to  all 
unknowns  in  each  event: 

37i/8XTi  =  [3R(S,Tiy8XTi  +  8R(Ti,Ri)/8XTi]/Vw 


=  [(Xji  -  Xs)/R(S.Ti)  +  (Xti  -  XRi)/R(Ti,Ri)lA^w  X->Y,Z 
aTi/8Xs  =  [aR(S,Ti)/aXs]/Vw  =  UXs  -  XTi)/R(S.Ti)]/Vw  X->Y,Z 


(A-4) 


8  71/8XRi  =  [8R(Ti,Ri)/8XRj]/Vw  =  [(Xrj  -  XTi)/R(Ti.Ri)]/Vw  X->Y,Z 

The  vector  L  contains  the  difference  between  observations  computed  with  the  model  and 
tl^  observed  data. 

After  the  normal  equations  are  formed,  a  priori  uncertainty  on  ship  positions  are  used  to 
constrain  the  solution.  Weights  are  added  to  the  normal  equations,  to  the  diagonals  of  N22, 
computed  from  ship  uncertainties. 
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The  coordinate  corrections  for  ship  positions  (S  and  Ri)  are  formally  eliminated 
algebraically  giving 

(Nil  -  Ni2(N22  +  P2)‘^N2i)Xi  +  (Ui  -  Ni2(N22  +  P2)'^U2)  =  0  (A-5) 

Each  event  j  is  processed  in  this  fashion  and  the  normal  equations  are  accumulated  for  the 

trackfc- 

NkXi  +  Uk  =  0 


where 

Nk  =  XjLi  1  ■  Ni2(N22  +  P2)‘^N2l)j 

Uk  =  X"=i  (Ul  -  Ni2(N22  +  P2)-lU2)j 

A  track  solution  for  the  coordinates  of  the  transponders  is  given  by: 

Xi=Nk-lUk 


(A-6) 


(A-7) 


Track  normal  equations  are  accumulated  and  a  flnal  solution  is  obtained  after  all  data  is  processed 
as 


xi=(Px+x”.i  (A-s) 

where  Px  is  a  weighting  matrix  developed  from  a  priori  uncertainty  on  the  location  of  the  array. 


ALTERNATIVE  OBSERVATION  MODEL 


The  observation  equation  (A-1)  presented  earlier  is  an  explicit  formulation  in  terms  of  time 
and  relies  on  a  knowledge  of  the  average  sound  speed  in  water  to  convert  metric  distances  into 
transit  times.  The  scale  of  the  resulting  network  is  directly  influenced  by  the  choice  of  this  velocity 
Vw- 


An  alternative  procedure  is  to  treat  the  observation  (transit  time)  as  an  implicit  function  of 
the  coordinates  for  the  transducer  and  transponders  in  the  following  way:  Consider  the  time  of 

transit  7i  associated  with  an  emitted  pulse.  The  angle  of  depression  asTj  at  the  emit  time  and  the 

angles  of  depression  aRjTi  at  the  receive  times  are  functions  of  the  positions  of  the  transducer  at 
points  S  and  Ri  and  the  positions  of  the  transponders;  i.e., 

aSTi  =  a(Xs,Ys,Zs,XTi,YTi,ZTi) 

and  i=l, 2,3,4  (A-9) 

ttRiTi  =  a(XRi,YRj,ZRi,XTi,YTi,ZTi) 
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where  Xs,Ys,Zs,XRj,YRi,ZRj,XTi,YTj,  and  Zjj  are  known  approximately. 

Using  a  water  velocity  profile  as  a  function  of  depth,  use  ray  tracing  to  determine  the  time 
of  transit  rr(S,Ti)  and  TrfRi.Ti)  using  the  angles  asXi  and  otRjTj  to  initiate  the  process.  Then  the 
observed  time  T\  is  related  to  Tx  by 

Ti  =  Tr(S.Ti)  +  TxiRiJi)  +  ATi  i  =  1,2,3,4  (A-10) 

The  equation  (A-10)  for  Ti  represents  the  (implicit)  observation  equation. 

In  developing  the  estimation  equations  (A-3)  using  this  approach,  the  partial  derivatives 
(A-4)  of  Ti  with  respect  to  any  of  the  coordinates  Xs,Ys,Zs,XRj,YRj,ZRj,XTi,YTi,ZTi  are  found 
by  a  variation  of  parameters  procedure.  That  procedure  is  based  on  determining  the  changes  in  the 
computed  observation  7)^  as  a  result  of  a  sequence  of  small  perturbations  in  the  initial  coordinate 
values.  These  are  then  used  to  develop  the  partial  derivatives  by  the  use  of  appropriate  ratios. 
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